Asymptotic near nucleus structure of the electron-interaction potential in local 

effective potential theories 
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In local effective potential theories of electronic structure, the electron correlations due to the Pauli 
exclusion principle, Coulomb repulsion, and correlation-kinetic effects, are all incorporated in the 
local electron-interaction potential v ee (r). In previous work, it has been shown that for spherically 
symmetric or sphericalized systems, the asymptotic near nucleus expansion of this potential is 
u ee (r) = « ee (0) + j3r + 0(r 2 ), with u ee (0) being finite. By assuming that the Schrodinger and 
local effective potential theory wave functions are analytic near the nucleus of atoms, we prove 
the following via Quantal density functional theory (Q-DFT): (i) correlations due to the Pauli 
principle and Coulomb correlations do not contribute to the linear structure; (ii) these Pauli and 
Coulomb correlations contribute quadratically; (iii) the linear structure is solely due to correlation- 
kinetic effects, the contributions of these effects being determined analytically. We also derive 
by application of adiabatic coupling constant perturbation theory via Q-DFT (iv) the asymptotic 
near nucleus expansion of the Hohenberg-Kohn-Sham theory exchange v x (r) and correlation v c (r) 
potentials. These functions also approach the nucleus linearly with the linear term of ^(r) being 
solely due to the lowest-order correlation kinetic effects, and the linear term of v c (r) being due solely 
to the higher-order correlation kinetic contributions. The above conclusions are equally valid for 
systems of arbitrary symmetry, provided spherical averages of the properties are employed. 

PACS numbers: 31.15.Ew, 71. 15. Mb, 71.10.-w 



I. INTRODUCTION 

The basic idea underlying all time-independent local 
effective potential theories of electronic structure such as 
Hohenberg-Kohn-Sham density functional theory (HKS- 
DFT) [l[ and Quantal density functional theory (Q-DFT) 
, is the mapping from a system of N electrons in some 
external field F ext (r) =-yo(r), where v(r) is the ex- 
ternal potential, to one of noninteracting Fermions with 
equivalent density p(r). For example, the potential en- 
ergy of the electrons in the external field of the nucleus 
of atoms is v(r) = —Z/r, where Z is the atomic number. 
The Schrodinger equation for the electrons is (in a.u.) 



= (T + V + L7)$ = EV, 



(1) 



T 



where 

2 Si^j ]7 



-sEiV?, v ^ = EMn), u = 

the wavefunction, and E the energy. 

The density of the electrons is p(r) = (\I f |/5|\f r ), where the 
Hermitian operator p — <5(rj — r). The corresponding 
Schrodinger equation for the model electrons (S system) 
is 



correlations are those due to the Pauli exclusion princi- 
ple and Coulomb repulsion, and the correlation-kinetic 
effects that represent the difference in kinetic energy be- 
tweeen the interacting and noninteracting systems. The 
highest occupied eigenvalue e m of Eq. ([2]) is the nega- 
tive of the removal energy whether it is the ionization 
potential of atoms and molecules or the work function of 
metals UMM- 

It is evident that any rigorous knowledge of the ex- 
act structure of the potential v ee (r) , and an understand- 
ing of how the different electron correlations contribute 
to this structure, is critical to the accurate construction 
of the model system. In our previous work [5] we had 
shown that v ee (r) is finite at the nucleus for spherically 
symmetric system. (For more general case of arbitrary 
symmetry, see Ref. Q.) We assume that the asymptotic 
near nucleus structure of v ee (r) for spherically symmetric 
systems is of the general form 



Vee(r) = v ee (0) + f3r + ^r 2 



(3) 



1 



V 2 +v(r)+v ee (r)] ( f> i (r) = e i <p i (Ty,i = l,...,N. (2) 



The S system wavefunction is the Slater determinant 
so that p(r) = = Ei$Wi(r). The 

v ee (r) is the local electron- interaction potential in which 
all the many-body correlations are incorporated. The 



In Eq. and the following, r is defined as the elec- 
tronic coordinate relative to one of the nuclei in the sys- 
tem. v ee (0) represents the value at electron- nucleus coa- 
lescence. In this paper we prove the following via Q-DFT: 
(i) There are no contributions to the linear term due to 
correlations arising from the Pauli exclusion principle and 
Coulomb repulsion; (ii) These Pauli and Coulomb corre- 
lations contribute to the quadratic term; (iii) The linear 
term arises solely due to correlation-kinetic effects, with 



2 



the coefficient [3 being 



3 P (0) 



[t(0)-i s (0)], 



(4) 



where t(r) and t s (r) are the kinetic energy densities of 
the interacting and noninteracting S systems, respec- 
tively; (iv) Finally, employing adiabatic coupling con- 
stant perturbation theory [7J, l8( , we derive via Q-DFT the 
asymptotic near nucleus structure of the HKS-DFT ex- 
change v x (r) and correlation v c (r) potentials. Although 
the above results are derived specifically for spherically 
symmetric systems, they are also valid for arbitrary sym- 
metry provided a spherical average of these properties are 
taken. (The results of this paper have also been obtained 
within the framework of HKS-DFT [§].) 

To put the present work in context, we note that there 
had been a long standing controversy about the struc- 
ture of «ee(r) at and near the nucleus of atoms. The 
traditional way of determining the exact structure of 
v ee (r) employs methods that assume knowledge of the 
true density p(r). Then one determines, via numerical 
methods, the v KK (r) that generates the orbitals leading 
to this density [HI . The densities employed are those ob- 
tained from correlated or configuration-interaction type 
wavef unctions. While some work [ll|, [T^, EH on the he- 
lium atom showed v ee (r) to be finite at the nucleus, work 
by others [HI, EB] showed it to be singular there. This 
singularity is introduced by the inaccuracies in the nu- 
merical calculations, for instance, those in the density 
arising from the limitations of basis sets employed to 
solve the differential equations for the wave functions or 
one-particle orbitals ■ Another key reason [|| for the 
singularity of v ee (r) at the nucleus is the lack of satis- 
faction of the electron-nucleus coalescence condition by 
the approximate wave function from which the density 
is obtained. For the determination of its structure for 
few electron molecular systems [TtJ such as H2 and LiH, 
the v ee (r) is assumed to be finite at each nucleus. Addi- 
tional work on light atoms, showed it to be either finite 
flM [l9L I20I ] or to diverge at the nucleus [2l| . Expressions 
for v ee (r) at a nucleus have also been derived, but once 
again they are based on the assumption that [22j it is fi- 
nite there. In various approximations within HKS-DFT 
[H llljp, [H, H3, the v ee (r) also diverges at the nu- 
cleus |2l|7^ With the knowledge that v ee (r) is finite at 
the nucleus, what remained was to understand whether 
the structure of v Ke .M approaches the nucleus linearly or 



quadratically h| 



111 . |12| . Il3t |28| . It is the answer to this 
question that we provide in this work. 

We begin in Sec. II by giving the key equations of time- 
independent Q-DFT. In Sec. Ill we provide the asymp- 
totic near nucleus expansion of the interacting and non- 
interacting system wavefunctions, and the corresponding 
expressions for the equivalent densities. The proof that 
Pauli and Coulomb correlations contribute quadratically 
and not linearly to the near nucleus expansion of v ee (r) 
is given in Sec. IV. In Sec. V we prove that the contribu- 
tion to the linear term of v ee (r) is solely a consequence 



of correlation-kinetic effects, and derive an expression for 
this contribution. Combining the results of the previous 
two sections, we present in Sec. VI the expression for the 
near nucleus structure of v ee (r). Finally, by employing 
coupling constant perturbation theory within the frame- 
work of Q-DFT, we derive in Sec. VII the near nucleus 
structure of the HKS-DFT exchange and correlation po- 
tentials. Concluding remarks are made in Sec. VIII. 



II. BRIEF INTRODUCTION TO Q-DFT 

In time-independent Q-DFT 0, the potential 

v ee (r) is the work done to move a model fermion in the 
force of a conservative field T e ^' (r): 



Vee(r) = - fV^VMl', 



(5) 



where T e *f{v) = £ ee (r) + Z tc (r), Here S ee (r) is the 
electron-interaction field representative of the Pauli and 
Coulomb correlations, and Z tc (r) the correlation-kinetic 
field representative of those effects. For spherically sym- 
metric systems, the fields £ e e(r) and Z tc (r) are each 
conservative. Thus, v ee (r) — W ee (r) + W tc {r), where 
the work done W ee (r), Wt c (r) are separately path- 
independent. The field £ ee (r) is obtained via Coulomb's 
law from its quantal source, the pair-correlation density 
g(r,r'), as 



£ ee (r) = / dv'g(r,r') 



r — r 



r'|3 ' 



(6) 



where g(r,r') = (*| J2ijtj ^( r i~ r )^( r j~ r ')\^) / p( r ) ■ Tne 
pair-correlation density g(r,r') may also be written as 
g(r,r') = p(r') + p xc (r,r'), where p xc (r,r'), its nonlo- 
cal component is the Fermi-Coulomb hole charge distri- 
bution. The pair-correlation density g s (r,r') of the S 
system may also be expressed in terms of its local and 
nonlocal components: g s (r,r') — p(r') + p x (r,r'), where 
the Fermi hole charge p x (r,r') — — |7s(r, r')| 2 /2 j o(r), 
with 7 s (r,r') = J2ia 4% ( r )^( r ') the Dirac density ma- 
trix. Defining the Coulomb hole charge p c (r, r') = 
g(r, r') — g s (r, r'), the field S ee (v) may then be written in 
terms of its components as 



£ ee (r) =£„{r)+£ x {r)+£ c (r), 



(7) 



where the quantal sources of the Hartree £ff(r), Pauli 
£ x (r), and Coulomb £ c (r) fields are the density p(r), 
Fermi hole p x (r,r r ), and Coulomb hole p c (r,r'), respec- 
tively. It follows that W ee (r) = W H {r) + W x {r) + W c (r), 
where the Hartree Wh{t), Pauli W x (r), and Coulomb 
W c (r) potentials are the work done in the respective 
fields. The field Z tc (r) is defined as 



1 



■[z a (r)-z(r)], 



(8) 



3 



where z(r) is defined in terms of the single-particle den- 
sity matrix j(r, r') as 



i (r)=2^% i (r)/9r j , 



(9) 



with 



^•(r) = -ld 2 /dndr< + d 2 /dr J dr' i \ 1 (r, r')\^ r , (10) 

the kinetic energy density tensor. The single particle 
density matrix is defined as 

7(r,r') = TVy^ / dx 2 . . . dx N ty*(rs, x 2 , . . . , xjv) 

*(r's,X2,...,xjv). (11) 

The z s (r) is similarly defined in terms of the Dirac den- 
sity matrix 7 s (r, r'). 



III. NEAR NUCLEUS BEHAVIOR OF THE 
WAVEFUNCTIONS AND THE DENSITY 

We first present the near nucleus behavior of the wave- 
functions for the interacting system and the noninteract- 
ing S system, and of the density. As is customary and 
the case in prior work [3(J HH , we assume that the wave 
function near the nucleus is smooth. In other words, 
that the electronic wave function of a bound state for an 
atom or molecule is analytic except at the nuclear cusps. 
This assumption requires the stringent condition that the 
Hamiltonian be analytic except where two or more par- 
ticles coincide, a condition that may not be mathemat- 
ically rigorous. However, in mathematical physics, the 
assumption of smoothness of the wave function is stan- 
dard practice, and we proceed by making it. Hence, we 
write the ground state many-body wave function as 

*(r,X) =*(0,X)+a(X)r + 6(X)r 2 + ... 



^ [o im (X)r + h m (X)r 2 ]Yi m (r 
v=-l 
2 

]T b 2m (X)r 2 Y 2m (f) + . . . , 



(12) 



m=-2 



for small r, where r — r/r, and Yi m (f) is the usual 
spherical harmonics. We use X for the coordinates 
s, T2-S2, ■ ■ ■ , tnSn- The coefficients in the wavefunction 
satisfy the following relations: 



a(X) + Ztf(0,X) = 0, 
2& lm (X) + Za lm (X) = 0. 



(13) 



The above relations originate from the cusp of the wave- 
function at the nucleus, and were derived previously in 



Ref. [3Gj, l3l| . The expansion of the wavefunction is ob- 
tained by making use of the asymptotic expansion for 
small r 

1 °° 1 A i 

R =EE^i^(W), d4) 

1 11 (=0 m=-l 'i 

for i = 2, . . . , N. The Schrodinger equation Eq. (|TJ) may 
then be rewritten as 



m— — 1 



^r 1 1 (X)*(r,X) = £*(r,X), 



where 



4 N 1 
3m (X) = — V —Y^ m {h), 

6 i=2 r > 



(15) 



(16) 



and HgZi (X) is the Hamiltonian for N— 1 electrons in an 
external potential of —(Z — l)/r. The first relation in Eq. 
(fl~3| follows readily by equating the coefficients of 
terms, and the second one by equating the coefficients of 
r°Yi m (f) terms. 

The density is defined as 



p(r) = N J dX|*(r,X)| 2 , 



(17) 



where J d~K denotes Y] s J e?xi . . . dxjy. Bear in mind that 
we are considering spherically symmetric or sphericalized 
systems. Employing Eq. (|12p in conjunction with Eq. 
(fT3")) and the relation 



Yi m (r)Yv m >{r)dr = S w S mm ,, (18) 

we obtain the density as 

p(r)= N J dx(V(0,X)| 2 (l-Zr) 2 
+ 2r 2 i?e[**(0,X)&(X)] 

+ r * E ^l«im(X)| 2 )+0(r 3 ). (19) 



Evidently 

p(0) =N J dX|*(0,X)| 2 . (20) 
One then can see that the density behaves as 
p(r) = (1 - Zr) 2 p(0) +r 2 2N J <JXfle[#*(0,X)&(X)] 

r 1 i 

+r 2 Nj rfX ]T ^|aim(X)| 2 + 0(r 3 ). (21) 
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Once again, with the assumption of smoothness, the 
orbital of the S system may be written as 



(f>nim(r) = r l [ A nl + B nl r + C n ir 2 
+D rd r 3 + ...}Y lm {f). 



(22) 



Similarly the coefficients in the orbital satisfy the cusp 
relations: 



B n o + ZA n Q — 0, 
2B nl + ZA nl = 0. 



(23) 



The density of the S system, which is supposed to be 
equivalent to that of the interacting system, is related to 
the orbitals as 



(24) 



nlra 



In the preceding summation over the one-particle states 
and the remaining ones in the paper, the spin summation 
is implicitly included. Substitution from Eq. (|2"2"]l leads 
to 



nlm 

+ \B nl \ 2 r 2 + 2Re{A* nl C nl )r 2 ]\Y lm (f)\ 2 . (25) 
Now employing Eq. (f2"3"]l . we obtain 

I 



p(r) = p(0)(l-Zr) 2 +r 2 — 
+3^L4 n i 



47T 

0(r 3 



2^J?e« C n0 ) 



(26) 



The fact that p{r) is a radial function, together with Eq. 
(fl"8|) . has also been employed to obtain Eq. (|26|) . 



IV. ELECTRON-INTERACTION FIELD £ es (r) 
NEAR THE NUCLEUS 

In this section, we show that the electron-interaction 
field £ ee (r) vanishes at the nucleus. This can be under- 
stood on physical grounds. The pair-correlation density, 
g(r, r'), the quantal source of the field £ ee (r), is in general 
not spherically symmetric about the electron position at 
r. However, for an electron at the nucleus, this charge 
distribution is spherically symmetric. Hence, the field 
due to this charge as obtained by Coulomb's law, must 
vanish at the nucleus. This result can also be proved 
mathematically. 

For spherically symmetric systems, the field £ ee (r) can 
only have a radial component, and can only depend upon 
the radial coordinate, i.e. £ ee (r) = r£ ee>r (r). Thus, by 
employing the identity 



f ■ y = d/dr, 



(27) 



one has, from Eq. ((SJ), 



£ee,r{r) = - J fa' g{v ,v') 



„ d 1 



dr |r — r' 



(28) 



According to Eq. (fT2|) . obviously g(r, r') is well behaved 
(not divergent) at r = 0. Therefore, 



£ee,r (?*) \r— 



1 

E Y lm (f)Y* m (f') 



m— — l 

o. 



(29) 



(For open shell atoms, the central field approximation 
means taking the spherical average of the field. Thus, in 
this approximation, the field is again radial and depends 
only on the radial coordinate.) 

Now since £ ee (r) is the negative gradient of the scalar 
potential function W ee (r), it is evident that there can 
be no linear term in this function. Thus, there is no 
contribution to the linear term of v ee (r) due to Pauli 
and Coulomb correlations. These correlations contribute 
quadratically. As an example, see Ref. [32| . 

Arguments similar to those above can also be made 
to prove that the Pauli £ x (r) and Coulomb £ c {r) fields 
defined as 



£x,c( r ) = / dr 'PxA r ' r ')-< 



„/ 3 : 



(30) 



vanish at the nucleus. (By the notation £ x c (r) we mean 
the Pauli ^(r) and Coulomb £ c (r) fields due to the Fermi 
p x (r,r') and Coulomb p c (r,r') holes, respectively.) Once 
again, for example, for the electron position at the nu- 
cleus, the Fermi and Coulomb holes are spherically sym- 
metric about the electron. Hence, the corresponding 
fields must vanish there. Hence, these individual cor- 
relations contribute to the quadratic structure of v ee (f) 
via the work W x (r) and W c (r) [12]. 

The classical Hartree field which is defined as 



S H (r) = / drV(r') 



r>\3 ' 



(31) 



corresponding to the Hartree potential, u_y(r) = 
J G?r'p(r')/|r — r'|, also vanishes at the nucleus because 
the density p(r) is spherically symmetric. 

In conclusion, the electron-interaction field £ ee (r), and 
its Hartree £ff(r), Pauli £ x (r), and Coulomb £ c (r) com- 
ponents all vanish at the nucleus. 



CORRELATION-KINETIC FIELD Z tc (r) 
NEAR THE NUCLEUS 



We shall first present the calculation for the kinetic 
fields of the interacting system and the S system sepa- 
rately in Subsec. A and Subsec. B. In Subsec. C, we 
present that for the correlation-kinetic field Z tc (r). 
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The summations over the dummy cartesian indices are 
assumed in this section. 



Correspondingly, the expression of Eq. (I2ip for the den- 
sity can be rewritten as 



A. Kinetic field z(r) of the interacting system 

In order to make our presentation simpler, we rewrite 
Eq. HHJ as 



tf (0, X)(l - Zr) + a(X) ■ r(l - Zr/2) 
+& j (X)r i r J , (32) 



*(r, X) 



where we have made use of Eq. (|13l) . The relations be- 
tween ot(X) and ai m (X), and /3jj(X) and &(X), &2m(X) 
are simply those between the cartesian tensors and spher- 
ical tensors. By substituting Eq. ([32]) into Eq. (jlip . we 
have the structure of the 7(1", r') for small r, r' as 

7 (r,r') = p(0)(l - Zr)(l - Zr') 

Zr){l- Zr'/2)A* -r' 
+(1- Zr')(l- Zr/2)A-r 
+E* j (l-Zr/2)(l-Zr'/2)r i r' j 
+[D* j (l-Zr) + F? ij r l ]r' i r' j 
+[D ij (l-Zr') + F Uj r' l ]r i r j , 



(33) 



where 



Ai=N / dXV(0,X)a*{X), 



E i:i = N J dXai(X)a*(X), 
Dij = N J dX*(0,X)/%(X), 



F iit = N / dXa t (X)(3*(X). 



(34) 



(35) 



(36) 



(37) 



It is necessary to mention that the preceding expression 
for 7(1:, r') is not limited to the spherically symmetric or 
sphericalized systems, but valid in general. For systems 
with central-reflection symmetry, Eq. (f3"3")) reduces to 

7(r,r') = p(0)(l - Zr)(l - Zr') 

+ ££(1-Zr/2)(1-Zr72)r i r$ 

+ - Zr)r# + Ai(l - Zr')nr r (38) 

We note that the second, third and fourth terms on the 
right hand side of the above expression are absent in the 
literature [3l|, HH. Also for completeness, we give the 
following relations: 

Da = 3N J dX**(0,X)6(X), 

E»= W J dX±-\a lm (X)\ 2 . (39) 



p(r) = p(0)(l - Zr) 2 + -{E u + 2D u )r 2 + 0{r 3 ). (40) 

The derivatives in calculating the kinetic energy den- 
sity tensor via Eq. (TTU|) are tedious but straightforward. 
After performing them, we arrive at 



Z 2 (p(0)+A.r)^ 



_ 1 

Z 
~Tr' 



iy(r) = 2 Re 

-(2 - Zr){r i A j + r 3 A t ) + E l3 {\ - Zr) 



+2(F tjl + F jU )n 



+ 0(r 2 ). 



(41) 



The corresponding kinetic energy density is 

t(r) = t u {r) = ^[Z 2 p(0) + E vl ] 
2Z 

-—riDu+E^+Oir 2 ) 



(42) 



The correctness of the expression for t(r) is indepen- 
dently confirmed via its definition as 

t(r) = ±N J dX V **(r,X). V *(r,X). (43) 

With the result for Uj(r) given in Eq. (|4ip. one can 
obtain from Eq. @ the kinetic field z(r) to be 

z(r) = 2Z 2 p(0)^ - \z{ZDu + E u )- + 0(r), (44) 
r z 6 r 



z(r) = 2Z 2 p{0)-^(l + 2Zr/3) 

~Z[2t(0) + 3Du]-+O(r), (45) 
3 r 

where we have made use of the following relation derived 
from Eq. (|32jl: 



Eu = 2t(0) - Z 2 p{0). 



(46) 



Equation (|46p can also be employed to rewrite the ex- 
pression for the density in Eq. (|40|) as 



p(r) = p( )|l-2Zr+— [Z 2 p(0) + t(0)]r 2 
+ ^r 2 D u +0{r 3 ). 
This expression for the density is employed later. 



(47) 
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B. Kinetic field z s (r) of the S system 



C. Correlation-kinetic field Z tc (r) 



The kinetic energy density tensor for the S system is 
expressed in terms of the orbitals as 



*»,«( r ) = \ Re Y^ 



d_ 

dri 



C/m( r ) 



_d_ 

dr 3 



0nim(r) 



(48) 



Substituting Eq. into the above expression and mak- 
ing use of Eq. (H23")) . one has, 



ST E \ A m\ 2 {$ij ~ ZrSij ~ Zr irj /r) 



8tt 

n 

- 2Z zZ ^Re(A* n0 C n o)r irj /r + 0(r 2 ).(49) 

n 

The kinetic energy density for the S system correspond- 
ingly is 

t s (r) = t Siii {r) 

n 

+ i(3-4Zr)^A| A „ 1 | 2 + 0(r 2 ). (50) 

n 

This expression for t s (r) is independently confirmed via 
its definition 



(51) 



(We note that the leading term i s (0) obtained from Eq. 
(jnnj) has been obtained previously [34j|.) 

The calculation for the kinetic field z s .i(r) = 
2^-i s ,ij(r) is now straightforward, which yields 



z s (r)= 2Z 2 p(0)^-12ZY / ^-MKoCno)- 

* — ' 47T r 



4Z£f \A nl f-+0{r\ 

— ' 47T r 



(52) 



or, by the use of Eq. (|50j) . 

z s (r) = 2Z 2 p(0)^{l + 2Zr/d)-4Z 1 - 



'M0) 



+ 0(r). (53) 



Also by the use of Eq. (|50|> . the expression of the density 
for the S system in Eq. (f2"6")l can be rewritten as 



p(r) = p(0)^l-2Zr + —{Z 2 p(0)+t s (0)}r 2 
+ 2r 2 J2^MKoCno) + 0(r 3 ). 

n 

This expression for the density is employed below. 



(54) 



We now substitute Eq. (J45J) and Eq. ([53]) into Eq. 
and arrive at 



4Z r 



(t(0)-t a (0)) + A 



p(r) r 



+ 0(r). (55) 



Now the density of the S system is equal to that of the 
interacting system. Comparing Eq. (|4T|) and Eq. (|54")) 
leads to 

i- £ i?e(^ C„o) = ^t(O) - t.(0) + D«]. (56) 
Making use of the above relation, one obtains, from Eq. 

2 tc (r) = ^-^UO) t(0)} V - + 0(r). (57) 

In contrast to the electron interaction field £ ee (r), the 
correlation-kinetic field -Z tc (r) is nonzero at the nucleus. 
For an example, see Ref. [351 ) - In passing, we mention 
Ref. [36| for a discussion of an approximate relation of 
the property J dr[t s (r) — t(r)] to p(0). 



VI. BEHAVIOR OF w ee (r) NEAR THE NUCLEUS 

Combining the results in Eqs. (f2"9")l and ([57]) . one finally 
has, from Eq. ([5]), 

4Z r 

V vee(r) = pr[t(0) - t.(0)]- + 0(r), (58) 

which also indicates the result shown in Eqs. ([3]) and ([4]). 
The classical Hartree potential due to the density, vu(r), 
is usually separated out from v ee (r) and the remaining 
part is referred as the HKS-DFT exchange-correlation 
potential v xc (r). Since it is evident from the structure of 
the Hartree field (r) that there is no linear term in the 
Hartree potential vh(t), Eqs. ([3]) and Q imply that 

v xc (r) = v xc (Q) + r-7^[*(0) - t.(0)]r + 0{r 2 ). (59) 
3p(0) 

Our calculations hence show that the linear term of 
v ee (r) (or equivalently v xc (r)) arises solely from the 
correlation-kinetic effects. 



VII. BEHAVIOR OF THE HKS EXCHANGE 
v x {r) AND CORRELATION v c (r) POTENTIALS 
NEAR THE NUCLEUS 

In HKS-DFT, the exchange-correlation potential 
v X c{y) is usually split into an exchange ^(r) and a cor- 
relation v c (r) component. These potentials have been 
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provided a rigorous interpretation in terms of the Pauli, 
Coulomb, and correlation-kinetic fields via Q-DFT @, [| 
as follows 



v x (r) = / [£ x (r')-Z t „ A (r')]-dl', 



(60) 

»c(r) = - fW)- Z to , 2 (r')-Z to ,3(r') -...]■ < 

(61) 

where Z tc l (r), and Z tc2 (r), etc., are the lowest order and 
higher order components, respectively, of the correlation- 
kinetic field Z tc (r) as obtained via the adiabatic cou- 
pling constant perturbation scheme @, H|- Thus, we 
see that HKS exchange potential v x (r) is representa- 
tive not only of Pauli correlation but also of the lowest- 
order correlation-kinetic effects. The correlation poten- 
tial v c (r) is representative of Coulomb correlations and 
higher-order correlation-kinetic effects. 

At this stage, it also becomes clear that the entire 
derivation of this paper can be reformulated within adi- 
abatic coupling constant perturbation scheme 0, Q • In 
such a scheme one has, instead of the result in Eq. ([55]) 
forZ tc (r), 

A7, r 

z *" (r) = Mty [ts{0) ~ <A(0)] ; + ° (r) ' (62) 

where A is the coupling constant. Define t x (r) as the 
first order correction to t s (r) in the adiabatic coupling 
constant perturbation scheme, and t c (r) = [t(r) — t s (r)] — 
t x (r), we have from Eq. 



(63) 



*..i(r) = 3^.(0); +0(r), 



2t.,a(r)+ 2t e , 3 (r) 



" Z i c (0)~ + O(r).(64) 



3p(0) v r 

Thus the work done in these fields is 

4Z 

W toil (r) = W te ,i(0) + 3^y*x(0)r + 0(r 2 ), (65) 

and 

Wt , 2 (r) + W tci3 (r) + • • • = W 4c>2 (0) + W tc , 3 (0) + . . . 

+ ^)^c(0)r + O(r 2 ). (66) 

Since we have shown in Sec. IV that £ x (r) and £ c (r) 
both vanish at the nucleus, we obtain the expression for 
v x (r) and v c {r) near the nucleus as 

AZ 

v x , c {r) = «», c (0) + T7^^^(°) r + °( r2 )- ( 6? ) 
3p(0) 

Thus, the v x (r) approaches the nucleus linearly, the term 
arising solely due to lowest-order correlation-kinetic ef- 
fects. The v c (r) also approaches the nucleus linearly, 



the contribution to this term arising from higher-order 
correlation-kinetic effects. 

Finally, it is worth mentioning that, as is well known, 
there is another definition of the kinetic ener gy d ensity 
of the interacting and noninteracting systems [371 ] : 



t(r)= -|jV / dXtt*(r,X)v 2 *U-.X) 



*«(r) 



nlrn 



Klm( T ) V 2 <AnZm(f). 



(68) 



Though t(r) and i s (r) yield respectively the same total 
kinetic energies for finite systems as t(r) and i s (r) defined 
in Eqs. (03J) and (JsTJ) , they differ from t(r) and t a (r): 



t(r) = t(r) + - V 2 P(r), 
i,(r) = *.(r) + i V 2 PW. 



(69) 



The second term on the right hand side of the preceding 
relations, which are singular at the nucleus are the same 
for the interacting and non-interacting systems. There- 
fore, all the results in this paper hold equivalently in the 
alternative definition. For instance, one also has, in par- 
allel to Eq. ([591 . 



v xc {r) = v xc (0) 



4Z 
37(0) 



[t(0)-t s (0)}r + O(r 2 ). (70) 



VIII. CONCLUSIONS 

In conclusion we note (i) the electron-interaction po- 
tential v ee (r) of the model electrons is finite at the nu- 
cleus. Correlations due to the Pauli exclusion princi- 
ple, Coulomb repulsion, and correlation-kinetic effects all 
contribute to this value. In this paper we have proved 
the following via Q-DFT: (ii) that the function v ee (r) 
approaches the nucleus linearly; (iii) there are no contri- 
butions to this linear term due to Pauli and Coulomb cor- 
relations; (iv) the Pauli and Coulomb correlations, how- 
ever, contribute to the quadratic term of this function; 
(v) the linear term in the structure of v ee (r) is solely due 
to correlation-kinetic effects. This contribution is deter- 
mined analytically. We have also obtained by application 
of adiabatic coupling constant perturbation theory via 
Q-DFT (vi) the near nucleus structure of the HKS-DFT 
exchange v x (r) and correlation v c (r) potentials. These 
potentials also approach the nucleus linearly: the linear 
term of v x (r) is solely due to the lowest order correla- 
tion kinetic effects; the linear term of v c (r) is solely due 
to higher-order correlation-kinetic contributions. We fur- 
ther note that the results in this work are valid in general 
whether the system is an atom, molecule or a solid. 

In previous work [38[, it has been shown that in 
the classically forbidden region of spherically symmet- 
ric systems, the contribution to the asymptotic struc- 
ture of v ee (r) due to Pauli correlations is of 0(l/r), 
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that of Coulomb correlations of 0(l/r 4 ), and those of 
correlation-kinetic effects of 0(l/r 5 ). Thus, together 
with the present results, the structure of v ee (r) in the 
asymptotic near nucleus and classically forbidden regions 
of such systems is now known analytically. (For the ana- 
lytic asymptotic structure of v ee (r) in the classically for- 
bidden region at the semi-infinite metal-vacuum inter- 
face, and the contribution of each type of correlation to 
this structure, we refer the reader to Ref. @, 

Finally, the result derived also has bearing on the con- 
struction of approximate HKS-DFT electron-interaction 
energy functionals of the density. In numerical codes 
employing such approximate functionals, the structure 
of the corresponding v ee (r) is obtained as the functional 
derivative. The derivative, being dependent on the den- 
sity, exhibits a sharp cusp at the nucleus just as 



does the density. However, calculations [LJ, l32| based 
on accurate wavefunctions show the slope of the function 
v ee (r) at and near the nucleus to be extremely small. 
Thus, these approximate energy functionals ought now 
be modified so that their functional derivative has the 
correct near nucleus asymptotic structure. 

We conclude by noting that the results of the paper 
can also be derived via the HKS-DFT 9]. However, the 
advantage of working within the framework of Q-DFT is 
that the near nucleus structure of v ee (r) can be under- 
stood as a function of the separate electron correlations 
that contribute to it. 

ZQ acknowledges support from the Chinese National 
Science Foundation under Grant No. 10474001. The 
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dation of CUNY. 
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